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Multipole correlations of t2g-orbital Hubbard model with spin-orbit coupling 
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We investigate the ground-state properties of a one-dimensional f2g-orbital Hubbard model including 
an atomic spin-orbit coupling by using numerical methods, such as Lanczos diagonalization and density- 
matrix renormalization group. As the spin-orbit coupling increases, we find a ground-state transition from 
a paramegnetic state to a ferromagnetic state. In the ferromagnetic state, since the spin-orbit coupling 
mixes spin and orbital states with complex number coefficients, an antiferro-orbital state with complex 
orbitals appears. According to the appearance of the complex orbital state, we observe an enhancement of 
r^u octupole correlations. 
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The competition and cooperation between spin and orbital 
degrees of freedom in strongly correlated electron systems 
manifest itself in the emergence of various types of spin- 
orbital ordered and quantum liquid phases. '"^^ In general, 
among competing interactions involving spin and orbital, the 
spin-orbit coupling is supposed to be weak in 3d transition- 
metal oxides such as cupurates and manganites, while as we 
move to 4d and 5d electrons, the spin-orbit coupling becomes 
strong and responsible for magnetic, transport, and optical 
properties. When the spin-orbit coupling is dominant, spin 
and orbital are not independent, but instead the total angular 
momentum gives a good description of the many-body state. 
In fact, it has been suggested that Sr2lr04, in which Ir''+ ions 
have five electrons in triply degenerate orbitals, exhibits 
a novel Mott-insulating state with an effective total angular 
momentum Jeff =1/2 due to a strong spin-orbit coupling. 
In the limit of strong spin-orbit coupling, the ground-state 
Kramers doublet at a local ion can be described by an isospin 
with Jcff=l/2.^'^* The exchange interaction among isospins 
can lead to a variety of ordering and fluctuation phenomena 
of spin-orbital entangled states. 

When we move to heavy-element /-electron systems such 
as rare-earth and actinide compounds, the spin-orbit coupling 
is large comparing with other energy scales. In such a case, 
we usually classify the complicated spin-orbital state from 
the viewpoint of multipole, which is described by the total 
angular momentum. Indeed, the multipole physics has been 
actively discussed in the field of heavy electrons.'*'^ A recent 
trend is to unveil exotic high-order multipole ordering. As an 
attempt to clarify multipole properties of /-electron systems 
from a microscopic viewpoint, we have numerically studied 
multipole correlations of an /-orbital Hubbard model on the 
basis of the j-j coupling scheme."' We believe that it is also 
important to clarify multipole properties in d-electron systems 
under the effect of the spin-orbit coupling. 

In this paper, we investigate multipole properties in the 
ground state of a one-dimensional t2g-orbital Hubbard model 
including the spin-orbit coupling by numerical methods. With 
increasing the spin-orbit coupling, the ground state changes 
from a paramagnetic state to a ferromagnetic state in terms of 
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the magnitude of the total spin. In the ferromagnetic phase, 
antiferro-dipole correlations develop even when the spin state 
is ferromagnetic due to the orbital contribution. On the other 
hand, the spin-orbit coupling induces a complex orbital state, 
in which real xy, yz, and zx orbitals are mixed with complex 
number coefficients. According to the complex orbital state, 
r4tj octupole correlations are enhanced. 

Let us consider triply degenerate orbitals on a one- 
dimensional chain along the x direction with five electrons 
per site. The one-dimensional f2g-orbital Hubbard model with 
the spin-orbit coupling is described by 

i,T.T' ,(7 i 

Pir^Piri + (^V^) PiraPir'a' 
+ iJm dljl,^,d,ra'd„'a 
+ (J72) dl,dl^,d.r'.'d.r'a, (D 

where dira (dlra) is annihilation (creation) operator for 
an electron with spin a (=t, i) in orbital r {=xy,yz, zx) at 
site i, and PiT(y=d\^„dira. The hopping amplitude is given by 
txy,xy=tzx,zx=t and zcro for other combinations of orbitals. 
Hereafter, t is taken as the energy unit. L; and Si represent 
orbital and spin angular momentum operators, respectively, 
and A is the spin-orbit coupling. U, U', J, and J' denote intra- 
orbital Coulomb, inter-orbital Coulomb, exchange, and pair- 
hopping interactions, respectively. We assume U=U'+J+J' 
due to the rotation symmetry in the local orbital space and 
J'=J due to the reality of the orbital function.'^' Throughout 
this paper, we set fcfcB=l. 

We investigate the ground-state properties of the model (1) 
by exploiting a finite-system density-matrix renormalization 
group (DMRG) method with open boundary conditions.'^' 
The number of states kept for each block is up to to=120, 
and the truncation error is estimated to be 10^"''^10^^. We 
remark that due to the three orbitals in one site, the number of 
bases for the single site is 64, and the size of the superblock 
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■^x^ Jy^ J z 


Tsg quadrupole 


0„=(1/2)(2J2-J2-J2) 
0,=(V3/2)(J2-J2) 


quadrupole 


*^yz—\V o/l)JyJz 

Ozx={V3/2)JzJx 






Oxy=(\/3/2)JxJy 


T2u octupole 


Txyz=(\/~^l(>)JxJyJz 


octupole 


T^=(m)(2Ji-j,j2-jij,) 

T^={V2){2Jl-JyJl-J^Jy) 
T^=(m(2Jl-JzJ^-Jpz) 


Fsu octupole 


Ti=iVT5mJxJ^-J^Jx) 

T^={VT5/6)iJy J2 - J2 Jy) 

T^=(VTE/6){JzJi-Jpz) 



Table 1. Definition of multipole operators up to rank 3. The overline on the 
product denotes the operation of taking all possible permutations in terms 



of cartesian components, e.g., JxJy=JxJy+JyJx 
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Fig. 1. Lanczos results for the four-site peiiodic chain, (a) The ground-state 
phase diagram in the (J, A) plane for U'=10. (b) The magnitude of the total 
spin in the whole system, (c) The correlation between spin and orbital in 
the single site, (d) The magnitude of the total angular momentum in the 
single site, (e) The charge density in each orbital. 



Hilbert space grows as m^x64^. To reduce the size of the 
Hilbert space, we usually decompose the Hilbert space into a 
block-diagonal form by using symmetries of the Hamiltonian. 
In the present case, however, the spin-orbit coupling breaks 
the spin SU(2) symmetry, so that we cannot utilize S'^^j as a 
good quantum number, where S^^^ is the z component of the 
total spin. Since we ignore eg orbitals among d orbitals, the 
total angular momentum is not a conserved quantity. The total 
number of electrons can be used as a good quantum number. 
Thus, since DMRG calculations consume much CPU times, 
we supplementally use a Lanczos diagonalization method for 
the analysis of a four-site periodic chain to accumulate results 
with relatively short CPU times. 

Let us first look at Lanczos results for the four-site periodic 
chain. In Fig. 1(a), we show the phase diagram in the ( J,A) 
plane for U'=10. The phase boundary is determined by the 
magnitude of the total spin S^^^. As shown in Fig. 1(b), Sj^t 
is almost zero for small A, indicating a spin-singlet ground 
state. As A increases, we find a transition to a ferromagnetic 
state with finite S^^^. Note that even in the limit of large A, 
Sjot does not approach the maximum value 2(2h-1)=6, since 
the spin-orbit coupling mixes spin up and down states and 
the complete ferromagnetic state is disturbed. In Fig. 1(c), we 
plot the correlation between spin and orbital in the local site 
Lioc • Sloe- At A=0, there is no correlation between spin and 
orbital. As A increases, the spin-orbital correlation develops 
and approaches one in the limit of large A, indicating totally 
parallel spin and orbital angular momenta. In Fig. 1(d), the 



magnitude of the total angular momentum in the single site 
Jfo^ is shown. At the transition point, 3f^^ exhibits a sudden 
increase, since the spin-orbit coupling stabilizes a large total 
angular momentum state at every local sites. Note again that 
Jj^Q^ does not reach the maximum value |(| + l)=x i'^ '■^^ 
limit of large A, since the total angular momentum is not a 
conserved quantity. Regarding the orbital state, we show the 
charge density in each orbital in Fig. 1(e). Due to the spatial 
anisotropy of orbitals, one hole is preferably accommodated 
in itinerant xy or zx orbitals in each site, while localized yz 
orbitals are doubly occupied. Measuring charge correlations, 
we find that holes occupy real xy or zx orbital alternately for 
small A (not shown). Namely, the ground state is a real orbital 
state. For large A, however, xy, yz, and zx orbitals are mixed 
with complex number coefficients by the spin-orbit coupling, 
leading to a complex orbital state. 

Now we move on to the analysis of multipole properties 
to clarify the ground-state properties from the viewpoint of 
multipole. We measure multipole correlation functions 



XrJq) 
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where XiY^ is a multipole operator with the symbol X of 
multipole for the irreducible representation in the cubic 
symmetry at site i. Here, we consider 15 types of multipoles 
including three dipoles (X=J), five quadrupoles {X=0), and 
seven octupoles (X=T), as listed in Table I.''** We evaluate the 
multipole correlation functions by DMRG calculations with 
chains of 16 sites. 

Figure 2 shows DMRG results of the multipole correlation 
functions at U'=10, J=2, and A=0 for the paramagnetic phase. 
Regarding dipoles, as shown in Fig. 2(a), the correlation 
has a peak at q=TT, which signals an antiferromagnetic state. 
Here, we notice that each of the dipole correlations exhibits a 
kink at q=Tr/2, while the kink corresponds to a peak for the 
Jy and Jz correlations. This kink structure originates in the 
spin-orbital SU(4) symmetry which realizes at a special point 
J=A=0.'^"'^* At the SU(4) symmetric point, correlations of 
spin Sj and orbital pseudospin Ti=^ Yl,r,r'.G dlra^TT'd-tr'a, 
where cr are Pauli matrices, coincide with each other and have 
a peak at q=n/2. With increasing J, the spin correlation of 
g=7r/2 grows and the peak of the spin correlation remains at 
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Fig. 2. DMRG results of muldpole correlations at U'=10, J=2, and A=0: 
(a) Tiu dipoles; (b) Fsg and T^g quadmpoles; (c) T4u octupoles; and (d) 
and F2u octupoles. 

q=Tr/2. On the other hand, the pseudospin correlation of q=TT 
is enhanced, and the peak position of the pseudospin correla- 
tion changes to q=TT. Note that for the orbital angular momen- 
tum hi, the correlation of the q=Tr component is enhanced as 
well. Thus, the orbital contribution to dipole leads to the peak 
of the Jx correlation at q=TT rather than g=7r/2. 

In Fig. 2(b), we find that for the 0„ and 0„ correlations, 
a sharp peak appears at q=0, since (Ou) and (0„) are turned 
out to be finite. We also find a peak at q=TT for the 0„, O^, 
and Oyz correlations, implying an antiferro-orbital state. For 
all quadrupoles, there occurs a kink at q=Tr/2 in similar to 
the case of dipoles, which is a trace of the SU(4) symmetry 
at J=A=0. As shown in Fig. 2(c) and 2(d), we also observe a 
kink at q=7r/2 for octupoles. 

In Fig. 3, we present the multipole correlation functions 
at U'=10, J=2, and A=0.5 for the ferromagnetic phase. At a 
glance, we find that the kink structure at q=TT/2 disappears 
for all multipoles. Concerning dipoles, as shown in Fig. 3(a), 
a peak appears at q=Tr. Namely, antiferro-dipole correlations 
become dominant even when the spin state is ferromagnetic 
due to the orbital contribution. In fact, the spin Si correlation 
has a peak at q=0, while the orbital Li correlation exhibits a 
peak at q=Tr (not shown). On the other hand, the quadrupole 
correlations are found to be almost flat, and we cannot see 
any fine structures signaling quadrupole ordering, as shown in 
Fig. 3(b). As for octupoles, we observe a significant enhance- 
ment of the octupole correlations of the q=TT component 
[see Figs. 2(c) and 3(c)]. The growth of the antiferro-octupole 
correlations reflects the stabilization of the antiferro-orbital 
state with complex orbitals. 

In summary, we have studied the ground-state properites of 
the i2g -orbital Hubbard model with the spin-orbit coupling 
from the viewpoint of multipole, by numerical techniques. 
The strong spin-orbit coupling induces a transition from the 
antiferromagnetic state to the ferromagnetic state. We have 
found that antiferro-dipole correlations develop even when 
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Fig. 3. DMRG results of multipole correlations at U'=10, J=2, and A=0.5: 
(a) Viu dipoles; (b) Fsg and F53 quadrupoles; (c) F4U octupoles; and (d) 
F511 and F2u octupoles. 

the spin state is ferromagnetic. Moreover, the complex orbital 
state appears, since the spin-orbit coupling yields the linear 
combinations of spin and orbital states with complex number 
coefficients. Accordingly, we observe an enhancement of the 
octupole correlations. It is an interesting issue to explore 
possible multipole ordering in Sd-electron Ir compounds with 
strong spin-orbit coupling. 
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